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Abstract 

In this paper, we calculated the next-to-leading order (NLO) new physics contributions to the 
mass splitting AMg d and the branching ratio B(B — ► X s ^) induced by the charged Higgs loop 
diagrams in the third type of two-Higgs-doublet models (model III) and draw the constraints on 
the free parameters of model III. For the model III under consideration, we found that (a) an 
upper limit \\tt\ < 1-7 is obtained from the precision data of AM Bd = 0.502 ± OMlps' 1 , while 
I At* I ~ 0.5 is favored phenomenologicaly; (b) for B — > X s j decay, the NLO QCD contributions 
tend to cancel the LO new physics contributions; (c) a light charged Higgs boson with a mass 
around or even less than 200 GeV is still allowed at NLO level by the measured branching ratio 
B(B -» X s j): numerically, 188 < M H < 215 GeV for (|A tt |, |A«,|) = (0.5,18); (d) the NLO 
QCD contributions tend to cancel the LO contributions effectively, the lower limit on Mjj is 
consequently decreased by about 200 GeV; (e) the allowed region of Mjj will be shifted toward 
heavy mass end for a non-zero relative phase 6 between the Yukawa couplings Xu and X^. The 
numerical results for the conventional model II are also presented for the sake of a comparison. 
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I. INTRODUCTION 



Among physical observables of B meson mixing and decays, the mass splitting AMg d 
and the branching ratio B(B — > Xsy) have been measured with high precision. Recent 
world averages as given in Refs. [H HI El are 

B{B -> X s7 ) = (3.34 ±0.38) x 10" 4 , (1) 
AM Bd = 0.502 ± 0.007 ps _1 , (2) 

which are in perfect agreement with the next-to-leading order standard model (SM) pre- 
dictions, for example, presented in Refs. 0, HI . Obviously, there is only a small room left 
for new physics effects beyond the SM. By comparing the theoretical predictions with the 
precision data, strong constraints on the parameter space of new physics models can be 
obtained. 

During the past decade, the B — > X s ^ decay has been studied in great detail in the 
SM and various new physics models. At present, the complete NLO calculations of this 
decay mode are available for the S M pi 0, H, ©, U|, for the conventional two-Higgs 
doublet models (2HDM) ©, Fl EJ for some supersymmetric models 

00, 

and 

for left-right symmetric models (lq]. The studies at leading order (LO) for model III 
0, 12, IE 3, f° r the top-two-Higgs-doublet model and the Technicolor models 0] 
are also available. From relevant theoretical calculations, we know that the rare decay 
B — > Xgj is very sensitive to the new physics contributions and has been used as the first 
test for new physics models. 

In Ref. 0, Borzumati and Greub studied the B —> X s ^ decay in NLO QCD in a class 
of models containing at least two Higgs doublets and with only one charged Higgs boson 
nondecoupled at low energy, and found the constraints on the mass of charged-Higgs boson 
for model II: the popular type-II two-Higgs-doublet model. In this paper, we will extend 
their work [ij to the case of model III ( the third type two-Higgs-doublet model0]), to 
constrain the charged Higgs boson mass Mh, as well as the the Yuakawa couplings A ti 
and Xbb after the inclusion of the NLO QCD corrections. 

The strength of the B® — B d mixing is described by the mass splitting AM# d . In the 
SM, AMfi d is strongly dominated by the box diagrams involving the heavy top quark and 
W gauge boson. In new physics models, the box diagrams where one or two W gauge 
bosons are replaced by other charged particles can also contribute to this quantity. The 
well measured AMg d therefore can be used to constrain the new physics models. 

In Ref.0, the authors presented a calculation of the mass splitting AM^ at NLO 
level in the conventional model II and found that the NLO corrections enhance the LO 
results by about 18%. In Ref. 0], the LO new physics contributions to AM^ have been 
calculated in model III and the constraint on the Yukawa coupling X u was also given by 



neglecting the large uncertainty of the non-perturbative parameter fs d y Bs d - We here 
will calculate the charged Higgs contribution to mass splitting AM^ at the NLO level 
in the model III, use the new high precision data to constrain \Vtd\ and \Xu\ and consider 

the effects of the large uncertainty of fg d J B Bd . 

This paper is organized as follows. In Sec. II, we describe the basic structures of the 
model III, give a brief review about the known constraints on model III presented in 
previous works. In Sec. Ill, we calculate the mass splitting AM^ at the NLO level in 
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model III, and draw the constraint on \X a \ by employing the new precision data. In 
Sec. IV, the NLO new physics contributions to the branching ratio B(B — > X s ^) in model 
III are calculated and analyzed in great detail. The numerical results for the conventional 
model II at NLO level are also presented as a comparison. The conclusions are included 
in the final section. 



II. STRUCTURE OF MODEL III AND CONSTRAINT 

The simplest extension of the SM is the so-called two-Higgs-doublet models [2^, 2^ . 
In such models, the tree level flavor changing neutral currents(FCNC's) are absent if 
one introduces an ad hoc discrete symmetry to constrain the 2HDM scalar potential and 
Yukawa Lagrangian. Let us consider a Yukawa Lagrangian of the formal 

A = VijQi,L$iU jt R + vFjQi^i D J,R + &i,L$zU jiR + ^Qu.o-'Dj.h + H.c, (3) 

where fa (i = 1, 2) are the two Higgs doublets, 0i j2 = iT 2 4>l 2 > Qi,L {U^r) with i = (1, 2, 3) 
are the left-handed isodoublet quarks (right-handed up-type quarks), Dj R are the right- 
handed isosinglet down-type quarks, while % '■ and ^ '■ = 1,2,3 are family index ) 
are generally the nondiagonal matrices of the Yukawa coupling. By imposing the discrete 
symmetry 



>2- 



A-> -Di,Ui^TUi (4) 



one obtains the so called Model I and Model II. In Model I the third and fourth term 
in eq.fJHJ) will be dropped by the discrete symmetry, therefore, both the up- and down- 
type quarks get mass from Yukawa couplings to the same Higgs doublet fa, while the fa 
has no Yukawa couplings to the quarks. For Model II, on the other hand, the first and 
fourth term in Eq.(|2I) will be dropped by imposing the discrete symmetry. Model II has, 
consequently, the up- and down-type quarks getting mass from Yukawa couplings to two 
different scalar doublets fa an d fa- 

During past years, models I and II have been studied extensively in literature at LO 
[ifi ] and NLO level 0, 0, 0, 0| and tested experimentally. The model II has been very 
popular since it is the building block of the minimal supersymmetric standard model. In 
this paper, we focus on the third type of 2HDM usually known as model III 17. 22l . 



In model III, no discrete symmetry is imposed and both up- and down-type quarks then 
may have diagonal and/or flavor changing couplings with <pi and 4>2- As described in [l7l |. 
one can choose a suitable basis (if , H 1 , H 2 , H^) to express two Higgs doublets [13] 



and take their vacuum expectation values as the form 



(5) 



<0i> = ( ,°/^ ), <02>=O, (6) 



v 



/V2 



where v = (v^2GV) = 2AQGeV. The transformation relation between {H° 1 H 1 ,H 7 



and the mass eigenstates (H°,h°,A°) can be found in |l7j. The H ± are the physical 
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charged Higgs boson, H° and h° are the physical CP-even neutral Higgs boson and the 
A is the physical CP-odd neutral Higgs boson. After the rotation of quark fields, the 
Yukawa Lagrangian of quarks are of the form |17| . 

£y 1 = ViA^iU^R + VijQi,L4>iDj t R + ^Qi,LhU hR + ^Q.joJ),,. + H.c, (7) 

where r^- D correspond to the diagonal mass matrices of up- and down-type quarks, while 
the neutral and charged flavor changing couplings will be 0] 1 



with 



^■neutral ~ S i (.charged ~ ? v CKMi (.charged ~ V CKMC, , 



* ~ V2M W Kv W 



where Vckm is the Cabibbo-Kobayashi-Maskawa mixing matrix 25[, i,j = (1,2,3) are 
the generation index. It is easy to see from Eq. Q that the the FCNC within the first two 
generations are naturally suppressed by the small quark masses, while a larger freedom in 
still allowed for the FCNC involving the top and bottom quark. The coupling constants 
Xij are free parameters to be determined by experiments, and they may also be complex. 

In the conventional model I and model II, the only additional contribution to the 
B — > X s 7 decay with respect to the SM comes from the charged Higgs boson-top quark 
penguin diagrams and depends on the mass of the charged Higgs boson, M H ±, and on 
tan/? = V2/V1, where v\ t 2 are the vacuum expectation values of (p\ t 2- From currently 
available studies at NLO level 0, one S e t to know the following main features 
of model II: 

• The charged Higgs penguins interfere constructively with their SM counterparts, and 
thus always enhance the branching ratio B(B —* X/y). The excellent agreement 
between the theory and experiments for the decay rate therefore leads to strong 
lower bound on the mass Mh±. One typical lower bound (99%C.L.) at NLO level 
as given in Ref.|9j is 

M H ± > 315GeV (10) 

for any value of tan /3. 

• The inclusion of NLO corrections shift the lower bound up by about 30% In 



other words, the lower bound on Mh± will become st rong er in the NLO level than 
that in the LO level. One lower bound as given in Ref. [12j is M H ± > 258 GeV using 
the LO calculation, but M H ± > 368 GeV using the NLO calculation. 

The lower bound on M#± from the measured branching ratio B(B — ► X^) depends 
very sensitively on small effects, and in particular on the way various errors are 
combined. The difference is usually about 100 GeV or even larger 12] . Since the 
theoretical error is significantly reduced at the NLO level, improving the calculation 
to the NLO has important effects on the lower bounds on Mh±. 



1 We make the same ansatz on the £^' D couplings as the Ref.0j. For more details about the definition 
of £ ' one can see Ref. 13. 
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• The bound on M H ± from B — » A, 7 decay is much stronger than those from other 
experiments. As shown in Fig. 4 of Ref. M, for example, the direct limit from LEP 
experiments is only M H ± > 78.6 GeV |2j, the ratio Rb = T(Z — > bb)/T(Z — > 
hadrons) is relevant only for very small tan f3, while rare B — > r decays constrain 
Mff± only for large tan/3. 

• The tan f3 dependence of the lower bound saturates for tan (3 > 5. 

For the model I, however, no bound on M#± can be obtained from B — ► A s 7 |27| . 
since the charged Higgs loops interfere destructively with the SM penguin diagrams and 
decouple for large tan/3. 

Althou gh m any phenomenolo gica l investigations have been done in the framework of 
model III |l6l 17, HH [l^, HHH^, the situation here is still not as clear as the model 
II, since there are much more free parameters in model III than in model II. As pointed 



in Ref. 17J, the data of K° — K° and B% — B® mixing processes put severe constraint 
on the flavor changing couplings involving the first generation of quarks. It is therefore 
reasonable to assume that the Yukawa couplings involving the u and d quarks are zero: 
Kj = ^dj = for j = 1, 2, 3. 

In Ref. [18], Chao et al., studied the decay b — > 57 at the leading order by assuming 
that only the couplings X a and A&& are non-zero. They found that the constraint on M H ± 
imposed by the CLEO data of b — > 57 can be greatly relaxed by considering the phase 
effects of Xu and A&&. The constraints from B° — B° mixing, the neutron electric dipole 
moment (NEDM), the Z°-pole parameter p and Rb were also considered in Ref.|l8|. The 



Chao-Cheung-Keung (CCK) scenario of model III |18( has following advantages: 

• Since one keeps only the couplings X tt and none zero, the neutral Higgs bosons 
do not contribute at tree level or one-loop level. The new contributions therefore 
come only from the charged Higgs loop diagrams with the heavy internal top quark. 

• The new operators O 9j i and all flipped chirality partners of operators O1...10 as 
defined in |2fiJ do not contribute to the decay b — > 57. 

• The free parameters in this model III are greatly reduced to Xu, Xbb and Mh±- 

In the following sections, we will calculate the NLO QCD contributions to the mass 
splitting AMs d and the branching ratio BiB — > X 8 j), to find the constraints on the 
parameters Xu, Xbb and Mh of model III. We will study the effects of the NLO QCD 
contributions in detail, and will also compare the results in model III with those in model 
II. 



III. B° - B° MIXING IN MODEL III 

B° — B° mixing is in general a FCNC process generated through weak interactions. At 
the lowest order of perturbation theory and in model III, the corresponding box diagrams 
which generates this process is shown in Fig. ^ The charged-Higgs boson contributions 
to B° — B° mixing at leading order were calculated long time ago [3(|. The NLO QCD 
corrections to B° — B° mixing was first presented in Ref. |3jJ for the SM, and in Ref. [5j for 
the conventional 2HDM: model I and model II. The possible constraints on model III from 
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the measured parameter Xd = AMb/Tb were studied, for example, in Refs. [l7L ll^. 13^ at 
the LO level. 



A. The basic formulae 

The strengths of the B® — B® mixing with q e (d, s) are described by the mass dif- 
ferences AMb = Mjj — M| where the subscripts H and L denote the heavy and light 
mass eigenstates respectively. The long distance contributions are estimated to be very 
small. The top box diagram is strongly dominant, while the charm and mixed top-charm 
contributions are entirely negligible. 

Recently, great progress have been made in experimental measurements. For AMb s , 
a lower limit of AMb s > 14. hps" 1 at 95%C.L. is available. For AMB d , however, it has 
been measured with high precision: the world average 0,0 is AM Bi = 0.502 ±0.007 and 
dominated by the results of B factories. At the end of the LEP-CDF-SLD era, AMb a has 
been measured with a relative precision of about 2.6% j33|. After including the B factory 
measurements, the precision is now 1.2% and high enough to constrain the new physics 
contributions effectively. 

On the theoretical side, the NLO theoretical prediction of AMs d is available now in 
the SM and in some new physics models beyond. In SM the mass difference AMs d can 
be written as [34| 

AM Bi = ^\m B \V td \ 2 (BB d f 2 B d )M^ V BS (x t ), (11) 



with 



Qtt 2 



S (x) = iX -] lx2 t/ -J^H*l (12) 



4(1 -x) 2 2(1 -x] 



B Bd = 5b» K 5 V)] 



-6/23 



3 

(5) 



QQ (aQ 

47T 



J: 



(13) 



Vb = [a s (/i 4 )] 6/23 x 



1 + a s (fit) f Si{x t ) + B _ j 
4tt \S (x t ) 1 



H In — t- + 7mn — ^ In 

2 Ml Jm0 dlnx t Ml 



(14) 



where x t = m t (fi t )/Ml 2 , 7 ^ = 4 and lm0 = 8 for SU(3) C , B t = 17/3 and J 5 = 
5165/3174 in the NDRscheme0, m B = 5.279 GeV is B° d meson mass, r} B = 0.55±0.01 
summarizes the NLO QCD corrections the function So(x t ) describes the dominant 

top-box contribution, fg is the B® meson decay constant, and B is the renormalization 
group invariant and non-perturbative parameter. There are a lot of works to estimate 
the values of fs d and Bb a in lattice QCD calculation and in QCD sum rule s 13311 . The 
definitions of the various quantities in Eq.([14|) can be found for example in Ref . |34j . Using 



m q (n) is the running q-quark mass in the modified minimal subtraction (MS) scheme at the renor- 
malization scale fi. For details, see Appendix 1X1 
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the input parameters as given in Appendix A, we find numerically that r] B = 0.553 and 
0.496 for fM t = 170 GeV or \i t = M w , respectively. The product rj B S (x t ) has, however, 
a very weak dependence on \i t at NLO level: the uncertainty is only 0.3% for 100 GeV 
<fit< 300 GeV. 

With a well measured AM Bd , one can determine \Vtd\ from Eq. ffTTj) in the framework 
of the SM. Using the input parameters as given in Appendix A, we find a limit on \Vtd\ 
from the measured AM Bd in Eq. (|TT| 



\V t 



td\ 



7.9 ± 1.6(f Bd y/ B Bd ) ± 0.2(m t ) ± 0.1( Vb ) 
(7.9 ± 1.7) x 10~ 3 . (15) 



x 10 -3 



where the uncertainties of f Bd y B Bd , m t and r/ B as listed in appendix A have been con- 
sidered and added in quadrature. In the following calculations, we will use this value of 
\Vtd\ as input. It is easy to see that the error of \V t d\ is almost completely determined by 



the uncertainty of the factor f Bd y B Bd . If this uncertainty can be decreased by a factor 
of two, we would find 

\V td \ = (7.9 ±1.1) x 10- 3 . (16) 



FigJ2] is the contour plot in f Bd \jB Bd — \Vtd\ plane obtained by the using the data 
AM Bd = 0.502 ± 0.007J9S" 1 . The shaded band in FigO shows the allowed region as given 



inEq.(fT5j). The solid line shows the f Bd \jB Bd dependence of \V t d\ and the width of the 
line shows the effect of uncertainties of all other quantities appeared in Eq. (fTT|) . 

B. Mass splitting AM B , in model II and III 

In the two Higgs doublet models, the charged Higgs boson contributes to the mass 
splitting AM Bd . In Ref . fl8j . the authors calculated the new physics contribution to AM Bd 
in model III at the leading order and presented the constraint on the X u — M B plane by 
using the measured Xd = AM Bd /T B . But they did not considered the effects of the large 



uncertainty of non-perturbative parameter f Bd y B Bd and the new physics contribution to 
the parameter rj B . 

Since AM Bd has been measured with very high precision, we use this quantity directly 
instead of the parameter Xd in our calculation. We will calculate the charged Higgs 
contribution to mass splitting AM Bd at the NLO level by extending the work in Ref.p 

to the case of model III. We will consider the effects of the large uncertainty of f Bd 
In the framework of 2HDM, the NLO mass difference AM Bd can be written as 

G 2 

AM Bd = ^m B M^\V td \ 2 (B B J B )r ]B (x t ,y t )S 2HDM (x t ,y t ), (17) 




67T 

with 

VB(x t ,y t ) = a s (M w 



,6/23 



1 as{M w ) ( D2HDAl{Xt,yt) _ ^ 
47T \S 2 HDM(x t ,yt) 
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and 



D2HDM 



[S {x t ) + S WH {x t ,y t ) + S HH (x t ,y t )} 
D S M(x t ) + D H (x t , yt), 



(19) 
(20) 



where x t = fn 2 (Mw)/M^ v and y t = m 2 (Mw)/M] I , and the high energy matching scale 
has been chosen as n = My/- The functions DsM(x t ) and Dn(xt,yt) in Eq. (j2T?j) describe 
the SM and new physics part of the NLO QCD corrections to the mass splitting AM^ 





D SM (xt) = C F [L^ SM \x t ) + 3S (x t )] + C A [L^ SM \x t ) + 5S (x t )] 
D H (x t ,y t ) = C F [L {1 ' H) (x t ,y t ) + 3 (S WH (x t ,y t ) + S H H(x t ,y t ))] 

+ C A [L^ H \x t ,y t ) + 5 (S WH (x t ,y t ) + S HH {x t ,y t ))} , 



(21) 



(22) 



where Cp = 4/3 and C A = 1/3 for SU(3)c- The function S (x t ) describes the dominant 
top-box contribution in the SM and has been given in Eq. (f*H?|) . The functions SwH(x t , yt) 
and SnH^Xt, yt) denote the new physics contributions from the box diagrams with one or 
two charged Higgs involved P 



SwH(x t ,y t ) = \Y\ 



Vt x t 



(2x t - 8y t ) Myt) 

(1 - yt) 2 (yt - x t ) 

-2x t 



+ 



6x t ln(x t ) 



1 - x t ) 2 (y t -x t ) 



yt)(i-x t ) 



S HH (y t ) = \Y\ 



yt x t 



l + y t 2y t \n[y t ] 



+ 



(23) 
(24) 



[l-ytf {l-ytf 

And finally, the function L^' H ^> {i = 1, 8) describe the charged Higgs contribution 

L^ H \x u y t ) =2\Y\ 2 WH^(x u yt) + 2\Y\ 2 Z>H®{x t , y t ) + \Y\ 4 HH^(y t ) . (25) 

The explicit expressions of complicated functions WH^(x t , yt), (x tl yt) and 

HH®(y t ) can be found in Ref.@. 

Following Ref.[4], we here use the symbols X and Y to denote the Yukawa couplings 
between the charged Higgs boson and quarks in the general 2HDMs. In the conventional 
model I and II, the couplings X and Y are real and given by 



X = -cot/3, K = cot/3 (Model I), 
X = tan/3, y = cot/3 (Model II). 



(26) 
(27) 



In the Model III where only the couplings \ tt and A^ are non-zero, the relation between 
the couplings (X, Y) and (X tt , Ay,) is simple 

X = -X bb , Y = X tt (Model III). (28) 

By using the input parameters as given in Appendix A, the SM prediction for AM# d 

is 



AM 



B d 



0.5061K ps- 1 



(29) 
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where the error comes from the uncertainty of parameter f Bd \J B Ba- 
in FigEl we show the M B dependence of AM^ in the model III , assuming \ u = 1. 
The region between two horizontal dot-dashed lines corresponds to the SM prediction 
as given in Eq.(|29|). The shaded horizontal band shows the world average AM Bd = 
0.502 ± 0.007ps _1 , and its width corresponds to the error. The short-dashed, solid and 



dashed curves in this figure show the model III predictions for f Bd y B Bd = 0.19, 0.23 and 
0.27, respectively. 

From the well measured physical observable AM Bd , one can find the constraint on 
| Ah | in model III. In FigEJ we show the |A#| dependence of AM Bd in the model III for 



f Bd y B Bd = 0.19 (the lower three curves) and f Bd y B Bd = 0.23 (the upper three curves) 
and for M H = 200 (solid curves), 250 (short-dashed curves) and 300 GeV (dashed curves), 
respectively. The shaded band and the region between two horizontal lines are the same 
as in FigCIl From this figure, an upper bound on \X a \ can be read off 

|A«| < 1.7. (30) 

This bound is complementary to the constraint obtained from FigO 

FigEl is the contour plot of the mass splitting AM Bd in the \Xu\ — M B plane, where 
the measured AM Bd and the relevant input parameters as given in appendix A have been 
used. The area A in FigJS] will be allowed by the measured AM Bd within 2a errors ( 
i.e., AM B , = 0.502 ± 0.014ps _1 ), if we do not consider the effect of the uncertainty of 



f Bd yB Bd in drawing this contour plot. This region corresponds to the allowed region 
in Fig.2 of Ref.jl^], but much narrow than that one because of the great progress of the 
experimental measurement of AM Bd . If we consider the effect of the large uncertainty of 



fB d y B Bd = 0.23 ± 0.04 in our calculations, the areas A plus B in the \X t t\ — M B plane 
will be allowed by the measured AM Bd within la error, while the areas A, B and C will 
be allowed by the measured AM Bd within 2a errors. 



Since the large uncertainty of f Bd y B Bd dominate the contour plot, one should take it 
into account in the effort to limit the free parameters in model III. As discussed previously, 



the uncertainty of \V t d\ are strongly correlated with the uncertainty of f Bd yB Bd , we 



therefore consider the uncertainty of f Bd y B Bd but use the central value \V t d\ = 0.0079 in 
our calculation. 

For the sake of experimental searches, one prefers a relatively light charged Higgs 
boson. From the contour plot FigEl the region of \X tt \ > 0.7 is disfavored if we expect 
existence of a light charged Higgs boson, while the parameter space of 

|A tt | » 0.5 and M H « 200GeV, (31) 

is certainly allowed by the measured mass splitting AM Bd . The allowed region in the 
contour plot will become narrow along with further improvement of the data and reduction 



of the large theoretical uncertainty of parameter f Bd y B Bd . 

Since the new physics contributions to AM Bd depend on \Y | 2 and \ Y | 4 only, the charged 
Higgs contributions in the model II and model III will be the the same if we use the same 
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value of the Yakawa coupling \Y\ as input. In model II, we have Y = 1/ tan (3. The upper 
limit on \X a \ as given in Eq. (j3U|) can be translated to a lower limit on tan/?, 

tan/3>0.6, (32) 

as can be seen directly from Fig® where the upper and lower three curves corresponds 

to fs d \J BB d = 0.19 and 0.23, respectively. The solid, short-dashed and dashed curves 
the model III prediction for Mh = 200, 250 and 300 GeV, respectively. The shaded band 
and the region between two horizontal lines are the same as in FigJ3J One can also see 
from this figure that the new physics contribution become negligible for tan/3 > 5. 



IV. THE DECAY B -► Y s7 IN MODEL III 

In the absence of new light degrees of freedom, the new physics contributions to the 
rare decay B — > Y s 7 will manifest itself through the new contributions to the Wilson 
coefficients of the same operators involved in the SM calculation, or the new operators 
absent in the SM, such as operators with different chirality. The excellent agreement 
between SM theory and experimental data leads to strong constraints on many new physics 
models beyond the SM. 

In this section, we calculate the branching ratio B(B — > X s 7) in model III. Here, the 
operator basis in the SM and model III under study is same. The NLO QCD corrections 
will be included in model III by extending the calculations in Ref. 0]. As a comparison, 
we also give the numerical results in the model II where it is necessary. 



A. Effective Hamiltonian and operator basis 

In the framework of the SM, if we only take into account operators up to dimension 6 
and put m s = 0, the effective Hamiltonian for b — > s-y at the scale \i reads [H 

AC 8 

H eff = --^V ts V^ C ^)°i^)- ( 33 ) 

* i=l 

The operator basis introduced by Chetyrkin, Misiak, and Miinz (CMM) 3 are given by 

Ox = (s L ^T a c L ) (c Ll »T a b L ) } (34) 
2 = (sl^cl) (c l j% l ), (35) 



3 There are two popular operator basis used in literature. The standard basis was defined for example 
in Ref-ji^]- The second one__is the CMM basis , where the fully anticommuting 75 in dimensional 
regularization are employed 
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3 = (s Ll ,b L )J2(irq), (36) 

g 

04 = (s L l»T a b L )J2(qrT a q), (37) 

q 

O, = (s L ^^ p b L )Y^(q^YYq), (38) 

<? 

O e = {s Ll ^ vlp T% L )Y^{^l v l P T a q), (39) 

<? 

7 = Y^m^is^bnjF^, (40) 

° 8 = lfe fr^™** ( 41 ) 

where T a (a = 1, . . . , 8) stands for SU(3) C generators, g s and e are the strong and electro- 
magnetic coupling constants, L,R = (1 =p 7s)/2 for the left and right-handed projection 
operators, Oi and 2 are current-current operators, 3 — 6 are the QCD penguin op- 
erators, and O-j and 0§ are electromagnetic and chromomagnetic penguin operators. In 
Eq. (J41j) . m&(/i) is the running b-quark mass in the modified minimal subtraction (MS) 
scheme at the renormalization scale \x (see appendix A for details). 

B. NLO Wilson coefficients at the scale nw and ^ 

To the first order in a s , the effective Wilson coefficients at the scale nw = 0(Mw) can 
be written as 0] 

cfim) = c?-W) + ^^^ ef W) • (42) 

The LO Wilson coefficients at the matching energy scale My/ takes the form 

Ct cf W) = l, (43) 

Cf ef W) = 0, (z = l,3,4,5,6), (44) 

C° 7 > eS (fi W ) = C° >SM (M W ) + \Y\ 2 Cl YY (M w ) + (XY*) C° tXY (M w ) , (45) 

C 8 °' ef W) = Cl SM {M w ) + \Y\ 2 Cl YY (M w ) + (XY*) Cl XY (M w ) , (46) 
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with 



°7,SM 



3x 3 t 



SM 



U 7,YY 



^7,XY — 10 



2x 2 : -8xf - 5x 2 + 7x t 
24(x t - l) 3 
-xf + 5x 2 + 2x t 

8(x t - l) 3 ' 
-8y 3 - 5y 2 + 7y t 
12(^-1)^"^ ' 72(^-l) 3 
-hy 2 t + % t - 3 + (6y t - 4) ln^ 



4(x t - l) 4 

4(x t - l) 4 
3% 3 - 2y 2 



In + - 
In x t + - 



fa " 1) 



/n0 



12(2/* - 1)' 



24(y t - l) 3 



a 



,XY 



Vt 

4 



-y 2 + 4y t -3-21ny f 

(w - 1) 3 



(47) 
(48) 
(49) 
(50) 
(51) 
(52) 



where x t = mf/M w , Vt = m t/Mx, and these leading order functions have no explicit nw 
dependence. 

The NLO Wilson coefficients at the matching scale nw hi model III can be written as 

1 



Ci' e (f^w) 



l.cff/ \ 
7 K^W, 

l,eff / 



15 + 6 In 



2 
3 



E + ^\nr^ + \Y\ 2 E H 
M w 



a 



C 8 ,e (vw) 



(z = 2,3,5,6), 

Cr^uifJ-w) + l^| 2 C] YY ([iw) + (XY*) C 7XY (^ W ) 
C^skAvw) + \Y\ 2 Cl )YY (n w ) + (XY*) C 8 xY (^ w ) 



l,SM 



where for i — 7, 8 the functions on the right-hand side of Eqs. ()56I57|) are 

CIyy(vw) = W iiYy + M itY] 
CIxy(^w) = W i:XY + M^ X) 



111 M 2 + J < 

In ifk + T- 



ln^# 

m 2 
/% 
In^l 



(53) 

(54) 

(55) 
(56) 
(57) 

(58) 
(59) 
(60) 



The explicit expressions of functions W{j, M^j and T it j (i = 7, 8 and j = SM, YY, XY ) 
can be found in Ref.|4|, and also listed in Appendix B for the convenience of the reader. 
The new physics contributions to the B — > X s ~f decay are described by the functions 

and (PEOJ). These 



C"f(n w ) with i = (7, 8) and j = (YY, XY) as defined in Eqs. 
eight functions depend on the unknown Mjj and the well measured Mw and mt only, 
and show clear decoupling behaviour when Mjj approaches infinity, as illustrated in Figd 
for the four functions C 7 ' YY (Mw) and C 7xY (M w ) . Other four functions C® YY (M W ) and 
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C s ' XY (Mw) have the very similar decoupling behaviour. Numerically, the LO new physics 
functions C®j(Mw) are always negative and have the same sign with their SM counter- 
parts, while the NLO new physics functions C}j(M w ) are always positive and have the 
opposite sign with their SM counterparts. 

Of course, the new physics contributions to B — > X s ^ also depend on the size and sign 
of the couplings X and Y, as can be seen easily from Eqs. (j49H52j) and ()59I60|) . For the 
conventional model II, we have 

' y|2 = s^ XY " = 1 ' (61) 

from the definition of X and Y as given in Eq. (j27|) . For model III under study, we have 

\Y\ 2 = \X tt \ 2 , XY* = -\\ tt \\\ hb \e ie } (62) 

where 9 — 9b — 9 t is the relative phase of the coupling X u and A&& in model III. The sign 
of the second term in above equation can be negative or positive depending on the choice 
of 9. In numerical calculations, we generally set 9 = 0°, unless otherwise specified, and 
will study the effects of a non-zero 9 in the end of this section. 

Since the heavy charged Higgs bosons have been integrated out at the scale fiw, the 
QCD running of the the Wilson coefficients Ci(fiw) down to the lower energy scale fib = 
0(mb) after including the new physics contributions is the same as in the SM. 

For complete NLO analysis of the radiative decay B — > X s 7, only the Wilson coefficient 
Cf{nb) should be known to NLO precision, 

CfW = Cr\^) + ^C^(^) , (63) 



with 



C° 7 ' cS ^ b ) = r/iC 7 ' cff (M + j| (r/^ (64) 

i=l 

^l,eff/ \ 39 x cff 8 / 37 39 \ „1 e ff/ v 

CV (W = V 23 Cj (nw) + g [V 23 -V 23 ) C 8 ' (nw) 

Z297664 is 7164416 m 256868 37 6698884 39 \ , eff/ . 

+ V 23 V 23 H V 2A V 2a Co' (Uw) 

V 14283 ' 357075 ' 14283 ' 357075 J 8 

I 37208 (J& n^i 



4761 



8 

+ Yl [ e *vCl' efi (nw) +fi + hv + hvCl' eS (fi W )] V ai , (65) 
i=i 

where the symbol r\ is defined as r\ = a s (fiw)/c(s(fib), and the "magic numbers" a^, hi, e^, 
fi, ki, and Zj in Eq. (j63jl are listed in Table |U 

The remaining coefficients are only needed to LO precision, 

8 

tf'^M = E M a % for J = 1, • • • , 6 , (66) 

i=l 

= V%C° 8 > eS (» w )+J2hiV ai , (67) 
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TABLE I: The "magic numbers" appeared in the calculations of the the Wilson coefficients 
Ci{n) in the rare decay b — ► 57. 



i 


1 


2 


3 


4 


5 


6 


7 


8 


Cli 


14 

23 


16 

23 


6 

23 


12 

23 


0.4086 


-0.4230 


-0.8994 


0.1456 


hu 








1 


-1 














h 2 i 








2 
3 


1 

3 














h 3 % 








2 

63 


1 

27 


-0.0659 


0.0595 


-0.0218 


0.0335 


ha 








1 

21 


1 

9 


0.0237 


-0.0173 


-0.01336 


-0.0136 


hu 








1 
126 


1 

108 


0.0094 


-0.01 


0.001 


-0.0017 


h%i 








1 

84 


1 

36 


0.0108 


0.0163 


0.0103 


0.0023 


e» 


4661194 
816831 


8516 
2217 








-1.9043 


-0.1008 


0.01216 


0.0183 


Si 


-17.3023 


8.5027 


4.5508 


0.7519 


2.0040 


0.7476 


-0.5358 


0.0914 




9.9372 


-7.4878 


1.2688 


-0.2925 


-2.2923 


-0.1461 


0.1239 


0.0812 


u 


0.5784 


-0.3921 


-0.1429 


0.0476 


-0.1275 


0.0317 


0.0078 


-0.0031 


hi 


2.2996 


-1.0880 


3 
7 


1 

14 


-0.6494 


-0.0380 


-0.0185 


-0.0057 


hi 


0.8623 











-0.9135 


0.0873 


-0.0571 


0.0209 


7i7 


208 
243 


416 
81 


176 
81 


152 
243 


6272 
81 


4624 
243 


32 
3 


32 
9 



where the "magic numbers " hi and hji are also listed in Table HI Following Ref . |4J , the 
small coefficients C® ,e (fib), ■ ■ ■ > C§' e (fib) will be neglected in numerical calculations. 

C. Branching ratio B(B — ► X s j) in the SM 

The branching ratio of the inclusive radiative decay B — > X s 7 can be written as 



B(B -> X sl ) LO = B SL 
at the LO level, and 

B(B -> X sl ) NLO = B SL 



v t * s v tb 



V, 



cb 



Kf(z) 



(68) 



v; s v tb 



v, 



cb 



6a f 



-Kf(z)n(z) 



D\ 2 + A + A] 



(69) 



at the NLO level, where Bsl = (10.64 ± 0.23)% is the measured semileptonic branching 
ratio of B meson 0, a em = 1/137.036 is the fine-structure constant 4 , z = m p ° le / 'm? h ole = 
0.29 ± 0.02 is the ratio of the quark pole mass 5 , f(z) and k(z) denote the phase space 



4 Based on the analysis in Ref. 35], it is more appropriate to use = 137.036 instead of a"* = 
130.3 ± 2.3 in the study of B — > X s -f decay. 

5 Another choice of z is z = rn c (/i) / rn^ ole = 0.22 ±0.04, but we do not consider this issue here. For more 
details about quark mass effects in B — > X s ^/ decay, see Ref.|9(. 
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factor and the QCD correction 

/(*) = 1 " ; 
k(z) = 1 



for the semileptonic B decay 



; z 2 + 8z e _ z 8 



24z 4 log(z) 



3tt 



7T 



+ 



riVP 



77? 



2 • 



(70) 
(71) 



where 8g[ denotes the nonperturbative correction to the semileptonic B meson decay, 



°SL 



2 2 



1-4- 



/(*) 



The term D in Eq. (p9*|l corresponds to the subprocesses b 

D = Cf(ji b ) + V(ji b ) , 



S7 



(72) 



(73) 



where the NLO Wilson coefficient Cj S (fib) have been given in Eq . (|63|) . and the function 
V(fib) is defined as 



An 



i=i 



1 ,u.,.ir ■ 



m 



16 „o, off/ 



(74) 



where the functions r,- 



. , 8) are the virtual correction functions (see Appendix D 



of Ref. 01), 7j°7 eff are the elements of the anomalous dimension matrix which govern the 
evolution of the Wilson coefficients from the matching scale fiw to lower scale /if,. The 



values of 7j°7 eff have been given in the last line of Table El for details see Ref. |34|. The LO 
Wilson coefficients in Eq.([74j) have been given in previous subsections. 

In Eq.fl69|). term A is the the correction coming from the bremsstrahlung process 
b -> j37[ 



A 



71 



Re-jr 



0,off 



(75) 



The coefficients /y have been defined and computed in Refs.jy, |3j]]. We here use the 
explicit expressions of those relevant as given in Appendix E of Ref. [H • 

In order to relate the quark decay rate to the actual hadronic process, the nonpertur- 
bative corrections obtained with the method of the heavy-quark effective theory (HQET) 
should be included. The term A NP in Eq.fjHHJ) and the term Sg[ in Eq. (|7T|) denote these 
nonperturbative corrections scale as l/mf and 1/m^ 



ml 



c° 7 ' eff (^ 



with 



2 s NP 
+ ^Re 
mi 



c-JVP 
5 1 



c° 7 ' eff (n b ) 



Ai 9. 

Y " 2 A2 ' 

9 ' 



(76) 

(77) 
(78) 
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where A2 = {m B » — m B )/4: = 0.12 GeV 2 extracted from B meson mass splitting, while 
the parameter Ai = 0.5 GeV 2 has large uncertainty but the overall Ai dependence largely 
cancels in the decay rate with the Ai term in 5$[ as can be seen from Eq. fJBTJlh 

Now we are ready to present numerical results of the branching ratios in the SM and 
2HDMs, specifically in model III. For the numerical evaluations, unless otherwise specified, 
we use the central value of input parameters given in Appendix [X] For the values of the 
matching scale and low energy scale, we always take nw = M w , and m b /2 < fi b < 2m b . 

Using Eq.fJHHl)? we find the SM prediction of the branching ratio B(B — > X s ^) 



3-52™Gu 6 ) ± 0.08(^)™ 

±0.13(a s ) ±0.05(m t ) ±0.04 
(3.52 ±0.28) x 10" 4 , 



m b 



v; s v tb 



v; 



cb 



x 10" 



(79) 



where the major sources of errors are shown explicitly, and the individual errors are added 
in quadrature. The central value 3.52 x 10 -4 is obtained by using the central values of 
input parameters listed in Appendix A and setting fi w = M w and \i b = m b . In this paper, 
we only consider the effects of the uncertainties of those six parameters as specified in 
Eq.flZi. 



D. Branching ratio B(B — > X s ^) in model III 

In the model III considered in this paper, the branching ratio has been parameterized 
in terms of three parameters X u , X bb and M#. From the limit on \X a \ obtained from the 
measured mass splitting AM Bd , the magnitude of the coupling X u should be smaller than 
1 if we require the charged Higgs boson is relatively light, as can be seen from FigJSJ In 
this section, we always set \X tt \ = 0.5, unless otherwise specified. 

We first check the common allowed regions in \X tt \ — \X bb \ plane for given M B = 250 
GeV. Fig. |S1 is the contour plot in the |A ti | — \X bb \ plane by using the LO and NLO 
theoretical predictions in model III and the measured branching ratio at the 2a level: 
2.58 x 10~ 4 < B(B — > Xg'y) < 4.1 x 10" 4 . In this contour plot, the first allowed region 
is the shaded area close to the X-axis where both \Xu\ and \X bb \ are all small, which is 
obtained when the NLO theoretical prediction is employed. The second allowed region 
is the area between two dashed curves obtained when the LO theoretical prediction is 
employed. The third region is the region between two solid curves allowed when the NLO 
theoretical prediction is employed. 

We choose one point of (X u , X bb ) in each of three allowed regions of Fig|H] as typical 
choices 

CaseA: (A«, X bb ) = (0.5, 1), 
CaseB: (A ti , X bb ) = (0.5, 12), 

CaseC: {X tt , X bb ) = (0.5, 22), (80) 

which will be used as input in the following numerical calculations. Here the limit on \X tt \ 
from the measured AM# d as studied in last section has been taken into account. 
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Since the new physics contribution to the decay B — > X s ^ is incorporated through its 
correction to the Wilson coefficients C^'l(fiw), we would like to check the size and sign of 
the new physics parts and their relative strength to the SM part, to show the theoretical 
features of the NLO contributions and to draw the constraint on Mh by compare the 
theoretical predictions with the data. 



1. Case A 

We firstly consider the Case A: (A tt , A^) = (0.5,1). Using the input parameters as 
given in Appendix A and assuming Mh = 250 GeV, we find the numerical results of Cy' 1 

C°' cS (M w ) = -0.1952 -0.0057 +0.0730 = -0.1280, (81) 

cf^{Mw) \y7c^{M w ) (XY^)C^(M w ) 

C^ eS (M w ) = -2.3712 +0.1013 -1.7239 = -3.9942, (82) 

C^M^) \Y\ 2 C$*{M w ) (XY^C^(M w ) 

at the matching scale fiw — M w , and 

Cf(m b ) = s -0.3137 t % +0.0082 j % +0.0507 j s -0.0142 / 

C 7'f M ( m b) ^C} SM {m b ) Cj'f p (m b ) /\C} NP (m b ) 
= -0.2690, (83) 

at the lower scale fib = m b , where ACjj(fi) = as ^/ C\'^ (//) with j = (SM, NP) denotes 
the NLO QCD corrections to the corresponding Wilson coefficients. It is easy to see that: 



• As shown in Eqs. (J81l82J) . the new physics contributions proportional to \Y\ 2 and 
XY* have the opposite sign and will cancel each other to some degree. The net 
new physics contribution is relatively small and has the opposite sign with its SM 
counterpart. 

• The net NLO new physics contribution has the same sign and comparable in size 
with its SM counterpart. 

• At the low energy scale //& = 0(m b ), the LO and NLO new physics contributions 
to CY S (m b ) has the opposite sign, they will cancel each other. The total Wilson 
coefficient C^ s (m b ) as given in Eq.(|83J) remains negative but changing from —0.3055 
to —0.2690 ( about a 12% decrease in magnitude) after the cancellation due to the 
new physics part. 

• The new physics corrections are generally small in magnitude, since the couplings 
X and Y are all relatively small for case A. 

Explicit calculations also show that the new physics corrections to V mb , A and A terms 
are induced through the modification to the Wilson coefficients Cj'£ s (fib)- But for case 
A, such corrections are small in size and therefore the theoretical prediction for branching 
ratio B(B — > X s ^) in this case agree well with the data within la error, as illustrated in 
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FigEl The region between two dot-dashed lines and the shaded part in FigEl shows the 
measured branching ratio with la and la error, respectively. The dashed line and solid 
curve shows the central value of the SM and model III (Case-A) prediction, respectively. 

Obviously, the theoretical prediction of the Case A agrees well with both the SM 
prediction and the data because both \ u and X bb are small in size. But this case is not 
interesting theoretically, since the new physics effect is too small to be separated from the 
SM contribution through experiments. 



2. Case B 

Now we turn to the case B: (X tt , \ bb ) = (0.5, 12). Her only leading order contributions 
in both SM and model III are taken into account, the branching ratio B(B — > X s ^y) 
completely determined by the Wilson coefficient C 7 ,e ^(u, b ) as shown in Eq. (j68|) . 

Using the input parameters as given in Appendix |X] and assuming Mh = 250 GeV, we 
find numerically that 

Cj' cS (M w ) = -0.1952 -0.0057 +0.8754 = +0.6745, (84) 

at the matching scale nw = Mw, and 

C 7 ' eS (m b ) = y -0.3137 / t -0.0042 / y +0.6581 j = +0.3402, (85) 

^0, eff / \ ^O.efT/ \ ^0, off / \ 

CV,sm(™6) C 7 , YY (m b ) C 7 \ XY (m b ) 

at the lower scale \L b = m b . It is easy too see from the numbers as given in Eqs. (j84l85|) 
that 

• The new physics contribution proportional to XY* term is much larger in size than 
the one proportional to X. 

• At both energy scales Mw and m b , the net new physics contributions to C 7 ' e (//) 
are always positive and much larger in magnitude than its SM counterpart. The 
total LO Wilson coefficient C° 7 ' eE (m b ) changed from -0.3137 in the SM to +0.3402 
in model III. 

From the B — ► X^ decay, only the magnitude of C-r{m b ) instead of its sign can be 
constrained by the relevant data. The semileptonic decay B — > K^l + l~ is sensitive to 
the sign of C 7 , but the precision of the measurement is still not high enough to determine 
the sign of Cj. 

In Fig. we show the Mh dependence of the branching ratio for the Case B. The 
dashed line shows the LO SM prediction BiB — > X s '-f) = 2.61 x 10~ 4 , while the solid 
curve shows the model III prediction for Case B. It is interesting to note that the LO 
SM prediction is marginally consistent with the data within 2a errors, but the NLO SM 
prediction agrees perfectly with the data. 

One can also find the limit on M H , 228 < M H < 264 GeV, from the Fig. ITU1 directly 
But one should know that the values of the lower and upper limits on M# will change 
along with the variation of \ tt and X bb . The point here is that a relatively light charged 
Higgs boson is still allowed in model III even at the leading order. 
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3. Case C 



= -0.1952 


+ 1.5992 = 


+1.4040, 


= -0.0972 


+ 1.4815 = 


+ 1.3843, 


= -0.0229 


- 0.3652 = 


-0.3881, 


= -0.0209 


- 0.2372 = 


-0.2580, 



For Case C, the LO and NLO new physics contributions to the Wilson coefficients 
Ci(Mw) and C$(Mw) are rather large. For given (Att,A&&) = (0.5,22) and assuming 
Mh = 250 GeV, we find numerically that 

C°' eS (M w ) = -0.1952 + 1.5992 = +1.4040, (86) 
C°' cS (M w ) = -0.0972 + 1.4815 = +1.3843, (87) 
a s (M w ) lee , 

— G 7 y M w) 

as ^ w ' Cl^{M w ) = -0.0209 - 0.2372 = -0.2580, (89) 

47T 

where the first and second term denote the SM and the new physics contributions, re- 
spectively. Clearly, the new physics contributions are always much larger than their SM 
counterparts in magnitude. 

After the inclusion of new physics contributions, the coefficients Cj ,eS (Mys) and 
Cg' e {My/) become large and positive. For NLO contributions, the new physics parts 
are around ten times larger than the corresponding SM parts. Among the new physics 
parts, the contribution proportional to term XY* is absolutely dominant over the one 
proportional to \Y\ 2 , for example, 

\\ tt \ 2 Cl YY {M w ) = -0.0057, -\\ tt \\\ bh \Cl XY {M w ) = 1.605 (90) 

for Case C and Mh = 250 GeV, since |A&fc| is now 40 times larger than |A tt |. 
Put the things together, we find numerically that 

Cf{M w ) = y -0.1952 / s -0.0411 / ^ +1.5992 y ^ -0.6548 / 

C°f M {M w ) ACl 3M {Mw) C° 7 f P (M w ) ACI NP {M W ) 
= +0.7081, (91) 

at the matching scale \iyr = M w , and 

Cf{m h ) = y -0.3137 / s +0.0082 / t +1.2024 j % -0.3320 / 

c 7,sli m b) ACj fllf (mi) C°'f p (m b ) AC] NP {m b ) 
= +0.5649, (92) 

at the lower scale fib = m b . From the numerical values in Eqs. (j91l92J) . one can see that 

• In the SM, the NLO QCD contribution AC} SM (M W ) has the same sign with its LO 
counterpart Cj'g M {My/) and is about 12% of Cj'g M {My/)- After the QCD evolution 
from Mw to low energy scale m&, the NLO QCD part AC^ SM {m b ) changed its 
sign and tend to cancel Cy\?^(m fe ). But the NLO QCD part is now only 2.6% of 
C?'g M {mb), and practically negligible for fi b w m b . 

• In model III and at the matching scale My/, the NLO QCD contribution 
ACj NP (M W ) has the opposite sign with its LO counterpart and is as large as 41% 
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of C?>° P (M W ). Unlike the situation in the SM, The NLO new physics contribution 
to C-j in model III will cancel its LO counterpart effectively, as shown explicitly in 
Eq.flHU). 

• In model III and at the low energy scale m b , the NLO QCD contribution ACj NP {mb) 

and Cj^ p (mb) still have the opposite sign, but the ratio of these two parts is lowered 
from 41% to 28%. The NLO and LO parts still cancel effectively. 

• Through the QCD running from My/ to mb, the Wilson coefficient C% SM changed 
its value from —0.2363 to —0.3055, increased by about 29% in magnitude. For the 
new physics part, Cj S NP is decreased by 8%. 

The numerical values in Eqs. (j91l92|) are obtained by assuming Mh = 250 GeV. For 
different values of Mh, the SM contributions remain unchanged, but the new physics 
part as well as their sum Cj S (rrib) will change greatly, as illustrated in Fig. ^2 The 
horizontal dots and dot-dashed line in Fig. HT1 corresponds to C^f M (mb) = —0.3137 and 
(mb) = 0.0082, respectively. The short-dashed and dashed curve shows the new 
physics contribution Cj'^ p {rrib) and AC^ NP (m b ), respectively. Finally the solid curve is 
the sum of those four parts. Since the dominant part Cj^ p {rrib) is decreased rapidly 
for increasing M H , the NLO part AC] NP {m b ) changed slowly and the SM parts remain 
unchanged, the Wilson coefficient Cj S (mb) is changing rapidly from large positive value 
to small negative value within the range of 200 < Mh < 800 GeV as shown by the solid 
curve in Fig. and will approach the value of its SM counterpart for heavier charged 
Higgs boson. 

One should note that when Cj S (fib) becomes small, other previously "small" NLO 
terms V(iib), A(fib) and A may play an important rule. The typical numerical values of 
those terms appeared in the curly brackets of Eq. (p)9"j) and their sum R, 

R = \~D\ 2 + A(m b ) + A, (93) 

are listed in Table |HJ In the SM, the relative strength of individual terms are 

Cf(m b ) : V(m b ) = 1 : 0.079, (94) 
R : \D\ 2 : A(m b ) : A = 1 : 0.98 : 0.030 : (-0.011) (95) 

The Wilson coefficient Cj S (rrib) clearly dominate the total contribution in the SM, while 
the radiative correction A(rrib) and the non-perturbative correction A play a minor rule 
since they are small in size and also cancel each other. 

In model III, however, the situation is very different because Cf{ m b ) and all other 
terms appeared in curly bracket of Eq. (J69j) will be changed by the inclusion of charged 
Higgs contributions through the modified Wilson coefficients C^ S (mb) and C^ s (mb), as 
can be seen from the numerical results listed in Table HP and the curves shown in Fig. IT2l 
The new features are the following 

• Along with the increase of Mh, Wilson coefficient Cj S (rrib) decrease rapidly, while 
V(mb) increase with less speed. Consequently, the " dominant" term \D\ 2 decreases 
more fast than Cj S (mb) does and approaches the minimum of \D\ 2 = 0.002 for 
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TABLE II: The numerical values of the NLO terms appeared in the curly brackets in Ea. (|69|) 
and the branching ratio B(B — ► X s j) (in unit 10 -4 ) in the SM and the model III of Case C, 
for typical values of mass Mh (in units of GeV). The term R is the summation of the terms 
\D\ 2 ,A(m b ) and A. 





Cf(m b ) V(m b ) \D\ 2 A(m b ) A R 


BTZ 


SM 


-0.3055 -0.024 - 0.015z 0.1089 0.0033 -0.0012 0.1110 


3.52 


200 
250 
300 


0.7559 -0.246 + 0.043i 0.2615 0.0138 -0.0513 0.2240 
0.5649 -0.208 + 0.036i 0.1283 0.0092 -0.0333 0.1042 
0.4210 -0.180 + 0.024i 0.0591 0.0064 -0.0222 0.0433 


7.10 
3.30 
1.37 


400 

Model III 600 
800 


0.2229 -0.140 + 0.020i 0.0073 0.0037 -0.0102 0.0008 
0.0105 -0.096 + 0.009i 0.0074 0.0023 -0.0022 0.0075 
-0.0951 -0.074 + 0.002* 0.0286 0.0022 -0.0002 0.0305 


0.03 
0.24 
0.97 


1000 
2000 
3000 


-0.1553 -0.061 - 0.002* 0.0467 0.0023 0.0003 0.0493 
-0.2580 -0.037 - O.OlOi 0.0873 0.0028 -0.0001 0.0900 
-0.2827 -0.031 - 0.012*' 0.0986 0.0031 -0.0006 0.1011 


1.56 
2.85 
3.21 



Mh ~ 480 GeV, which leads to a negative R and consequently a negative branching 
ratio: an unphysical result. For more heavier charged Higgs boson, the new physics 
contributions becomes smaller and smaller, while the summation of individual terms 
in the curly brackets of Eq. (JrJ§j) restores to its SM value slowly, as illustrated by the 
solid curve in Fig. For Mh > 600 GeV, the terms A(m b ) and A remain basically 
unchanged. As shown in Fig|7J the new physics contributions approaches zero when 
M H approaches infinity, this is the so-called decoupling behavior of Higgs boson. 

• From Table HP and FigIT2*l one can infer that there are two regions of Mh allowed 
by the measured BiB — > X s j). This point can be seen more directly from Fig. 

In Fig fEfl we draw the Mh dependence of the branching ratio BiB — > X s 7) in model 
III of Case C: (X u , X b b) = (0.5,22). In both (a) and (b), the band between two dot- 
dashed horizontal lines shows the data within 2a errors. The shaded band shows the SM 
prediction and the error as given in Eq. (j79|) . The dash and solid curve corresponds to 
the LO and NLO model III predictions, respectively. In order to show the effects of NLO 
corrections, we here use the same values of (X a , X bb ) = (0.5, 22) as input for both LO and 
NLO theoretical predictions. The FigfEIb is a magnification of the light Higgs part of 
FigfLlk For given X u = 0.5, X bb = 22 and M H = 250 GeV, we find numerically 



B(B 



X si) NLO 



3-48±8:S|(/z 6 ) ± 0.08(£W +0 ' 42 



,m. 



-0.36 < 



+0.52 



0.58V a s. 



±0.26(mt) ±0.04 



v t * s v tb 



V, 



cb 



x 10" 



(3.48 ± 0.74) x 10 



-4 



(96) 



where the central value of branching ratio is obtained by using the central values of input 
parameters as given in appendix while the six major errors from the uncertainties of 
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those input parameters are added in quadrature. For asymmetric errors we use the larger 
value in making quadrature. 

For Case C, the allowed regions of Mh can be read off from FigIT3*l directly. 

432 <M H < 478 GeV (97) 

at the LO level, and 

236 <M H < 266 GeV, and M H > 1640GeV (98) 

at the NLO level. The first allowed region, a region interested by many physicists from 
the point of experimental searches, is shifted to the lower part by about 200 GeV because 
of the inclusion of NLO contributions. 

If we consider the effect of the theoretical error as given in Eq. (j96|) . the limits at the 
NLO level will become 

226 < M H < 285 GeV, and M H > 1120GeV. (99) 

The first allowed region has a weak dependence on the theoretical error, but the lower 
limit of the second allowed region of Mh is very sensitive to the theoretical error. 

FigHH is the contour plot in \Xbb\ — Mh plane obtained by using the LO and NLO 
model III predictions and the measured decay rate at 2a level, while assuming |A tt | = 0.5. 
The regions between two dashed curves and two solid curves are allowed when the LO 
and NLO theoretical predictions are employed. One can see from Figs ED and El that a 
relatively light charged Higgs boson in model III, say around or even less than 200 GeV, 
is still allowed by the measured branching ratio B(B — ► X s j). For example, the region 

188 <M H < 215GeV, (100) 

is allowed by the data if we set (\\tt\, I--WI) — (0-5, 18) in model III. This is a good news 
for future experimental searches! 



E. Comparison of the results in model II and III 

Between the conventional model II and the model III studied here, there is a direct 
transformation. The results of model III can be reduced to the results of model II by the 
substitution 

\ bb = - X -)• - tan (3 and A« = Y -> l/tan/3. (101) 

As a comparison, we here also present the LO and NLO model II prediction for the 
branching ratio BiB — > X s 7) and show the lower limits on Mh obtained from the data. 
For more details about the rare decay B — > X s ^ in the model II. one can sec for example 
Ref. ji| and references therein. 

By using the central values of input parameters as given in Appendix El and assuming 
tan f3 = 4, we find the numerical values of the NLO terms appeared in the curly brackets 
in Eq. (j69|) and the branching ratio B(B — > Xg'j) in the model II for typical values of Mh, 
as listed in Table UTT1 In contrast to the model III of Case C, the differences between the 
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TABLE III: The numerical values of the NLO terms appeared in the curly brackets in Ea, (|69|) 
and the branching ratio B(B — ► X^) (in unit 10~ 4 ) in the model II with tan (3 = 4, for typical 
values of mass M# (in units of GeV). The term R is the summation of the terms \D\ 2 , A(rrib) 
and A. 





Cf{m b ) V(m b ) \D\ 2 A(m b ) A R 


BK 


Ql\/T 


n n n9/i n m n i dsq n nn , ?' : ! n nm 9 n 1 1 m 

— u.ouoo — — U.U1D? u.iuoy U.UUOO — U.UU1Z U.111U 


1 ^9 
O.OZ 


200 
250 
300 


-0.4034 -0.004- 0.020?' 0.1661 0.0042 -0.0037 0.1666 
-0.3857 -0.007 - 0.019? 0.1547 0.0040 -0.0032 0.1555 
-0.3723 -0.010 - 0.019?' 0.1464 0.0039 -0.0028 0.1475 


5.28 
4.93 
4.68 


400 

Model II 600 
800 


-0.3540 -0.014 - 0.018i 0.1354 0.0038 -0.0024 0.1368 
-0.3345 -0.018 - 0.017?' 0.1242 0.0036 -0.0019 0.1259 
-0.3247 -0.020 - 0.016? 0.1188 0.0035 -0.0017 0.1207 


4.34 
3.99 
3.83 


1000 
2000 
3000 


-0.3192 -0.021 - 0.016?' 0.1158 0.0035 -0.0015 0.1178 
-0.3098 -0.023 - 0.015?' 0.1109 0.0034 -0.0013 0.1131 
-0.3075 -0.024- 0.015? 0.1098 0.0034 -0.0012 0.1120 


3.74 
3.59 
3.55 



values of the NLO terms in the SM and model II are small and approaches zero when M# 
approaches infinity. It is easy to see from Table ILTTl that the new physics correction to the 
branching ratio becomes less than 6% and 2% for M# = 1000 and 2000 GeV, respectively. 

If we consider the effects of uncertainties of input parameters, we find the branching 
ratios at LO and NLO level 



li 

LO 



3.61^(^)±0.08(B Si )l°;f 7 (^) 



±0.09(a s ) ±0.08(m t ) ±0.04 
(3.61 ±0.82) x 10~ 4 , 







')] 


( 


v cb 





x 10" 



(102) 



and 



B(B - X sl )» LO 



4.i3i^(^ 6 ) ± om(B SL )l:» 

±0.09(a s ) ±0.08(m t ) ±0.04 
(4.13 ±0.34) x 10~ 4 , 



+0.24 / m C' 



2(1 \ 





v t * s v tb 


')] 


I 


Vcb 





x 10" 



(103) 



where the central value is obtained for M# = 500 GeV, and the errors connected with 
the uncertainties of six input parameters are added in quadrature. 

Fig. ITB1 shows the M# dependence of the branching ratio B{B — > X^) in model II at 
LO and NLO level, assuming tan/3 = 4. The shaded band and the band between two 
horizontal dot-dashed lines shows the measured branching ratio B(B — > X s ^j) within la 
and 2a errors, respectively. The dashed and solid curve shows the LO and NLO model II 
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prediction. By comparing the theoretical predictions and the measured branching ratios 
at 2cr level, it is easy to read off the lower limit on Mh directly from Fig. 



M H > 357GeV 



(104) 



at LO level, and 



M H > 520GeV 



(105) 



at NLO level, if the uncertainties of the input parameters are not taken into account. If 
we consider the combined uncertainties as given in Eqs. ()102ll03|) the lower limit will be 

changed to 



at the LO and NLO level, respectively. 

Fig. El is the contour plot in tan (3 — Mh plane obtained by finding the minimum of 
the NLO model II prediction of the branching ratio BiB — > X 3 j), when varying the input 
parameters within their errors and the lower scale fib in the range of 2.4 < fit < 9.6 GeV, 
but fixing fiw = M\y- The excluded region is below the corresponding curves when the 
measured branching ratio at la (dashed curve) and 2a level (solid curve) is employed. 
For B(B -> X s j) < 4.1 x 10" 4 and tan/3 = 1, 10,20, the lower limit on M H is 315, 278, 
and 277 GeV, respectively. Because of the flatness of the curves shown in Fig. Hoi towards 
the higher end of Mh, the lower limits on Mh are very sensitive to the ways to deal with 
theoretical errors or the details of the calculations. The lower limit on Mh as illustrated 
in Fig^Jis consistent with that given in Ref. |4|. 

The major differences between model II and the model III studied here and the causes 
inducing such differences are the following 

• In model II, the region of tan/3 < 1 is strongly disfavored by the measured AM^ 
and other experiments. In the region tan/5 > 4, we have \Y\ 2 = 1/ tan 2 /3 <C 1 
and XY* = 1, the possible new physics contributions to the Wilson coefficients Cj t % 
are small in magnitude and have the same sign with their SM counterparts, and 
therefore strongly constrained by the excellent agreement between the data and the 
SM prediction for B — ► X s ^ decay. 

• In the model III, the new physics contributions proportional to \\tt\ 2 and |A#||A&&| 
have the opposite sign for 6 = 0° and will cancel each other. Even if |Att| should 
be smaller than 1 due to the strong constraint from the measured mass difference 
AMb,, as discussed in last section, the inclusion of new physics contribution still 
can change the sign of the Wilson coefficient C| ' (ra&) from negative in the SM to 
positive in model III. On the experimental side, current date still can not exclude 
the possibility of a positive Wilson coefficient Cj S (rrib). 

• Other NLO terms in the curly bracket of Eq. (j69|) also receive new physics correc- 
tions, and therefore the cancellation between individual terms in model III is very 
different from the pattern in model II. 

• The inclusion of NLO contributions will decrease the lower limit on Mh in model 
III. In model II, in contrast, the lower limit will goes up by including the NLO 
corrections. 



Mh > 298 or 350 GeV 



(106) 
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• In the model III studied here, a light charged Higgs boson with a mass around or 
even less than 200 GeV is still allowed at NLO level. In model II, however, such 
light charged Higgs boson seems impossible. 

F. Effects of a non-zero phase 9 

For the mass splitting AMs d , the new physics contributions in model III depend on 
|A tt | 2 and |A t4 | 4 , and therefore are independent of the phase of coupling X tt . 

For the B — > X s 7 decay, however, the dominant new physics contributions ( i.e. the 
third terms in Wilson coefficients C° 7 f{M w ) and C l 7 f{M w ) in Eas. lHoTO and (|5EI57|) ) 
depend on the relative phase 9 between X tt and as can be seen from Eq.(|B2|). In 
previous calculations, we always assume 9 — 0°. Here we will consider the effect of a 
non-zero phase 9. 

FigEI shows the Mh dependence of the branching ratio B(B — ► X s ^) in model III 
for case C (i.e. (|A t t|, |A#|) = (0.5,22)) and for 9 = 0° (solid curve), 30° (dots curve), 
60° (short-dashed curve) and 90° (dashed curve), respectively. The band between two 
horizontal dot-dashed lines shows the measured branching ratio within 2a errors. The 
shaded area shows the SM prediction with the error as given in Eq. (J79|) . It is easy to see 
from FigEI that the non-zero phase will strengthen the constraint on the mass Mh in 
case C. For 9 = 30°, the allowed region is shifted to 

271 < Mjj < 317 GeV. (107) 

For 9 = 60°, the lower limit on Mh goes up to 417 GeV. For 9 > 75°, the lower limit on 
M H is higher than 600 GeV. 

In Fig ll8l we show the 9 dependence of the branching ratio B(B — ► X s ^) in model 
III for the case A (the short-dashed curve) and case C (the solid curve) discussed in 
previous subsections, and for given Mh = 250 GeV. Obviously, only two narrow regions 
of 9, 0° < 9 < 20° and 340° < 9 < 360°, are allowed by the data for case C. For case A, 
however, almost the whole ranges of 9 are still allowed by the data of B — > X s ~f decay. 

V. SUMMARY 

In this paper, we calculated the new physics contributions to the mass splitting AM^ 
and the branching ratio B(B — > X^) induced by the charged Higgs loop diagrams in 
model III, and found the constraints on the parameters of model III by comparing the 
theoretical predictions with the high precision data. We focus on the effects of the NLO 
QCD corrections to these two physical observables. 

In section |Hj we gave a brief review about the structure of the general two-Higgs- 
doublet models and the phenomenological studies about such models. Following previous 
works [l8j], we assume that only the diagonal Yukawa couplings \ a and A^ are non-zero 
for the model III under consideration in this paper and study the effects of NLO new 
physics contributions. 

In section Hill we calculated the new physics contributions to the mass splitting AM Bd 
in model III at the NLO level. The magnitude of Yukawa coupling X tt is strongly con- 
strained by the precision data AM^ = 0.502 ± 0.007ps~ 1 . As shown in FigsEJandEl the 
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upper limit on X u is \X t t\ < 1.7, while the choice of \X a \ ~ 0.5 is favored by the measured 
AM Bd if one requires the charged Higgs boson to be light, say around 200 GeV. 

In section ITVl we calculated the new physics contributions to the rare decay B — > X s ^ in 
the model III. The NLO QCD corrections are taken into account here. The new physics 
contributions to the Wilson coefficients Cy' 8 (//) and the interference between the new 
physics parts and their SM counterparts are investigated. The new physics contributions 
to B — > decay and possible constraints on \\bb\ and Mh are calculated and analyzed, 
three typical cases for the choice of Yukawa couplings \X tt \ and \\^\ are studied in great 
detail. The common features and the differences between the conventional model II and 
the model III, as well as the effects of a non-zero relative phase 9 are also considered. 
From the numerical results, we found that 

• In the model III studied here, a light charged Higgs boson with a mass around or 
even less than 200 GeV is still allowed at NLO level by the measured branching 
ratio BiB — > X s ^) within 2a errors, as can be seen in Eq. (jl(JU|) and Figs JT31 and ITU 
In model II, however, such light charged Higgs boson seems impossible. 

• The inclusion of the NLO QCD contributions will decrease the lower limit on M# 
by about 200 GeV in model III. In model II, in contrast, the lower limit on Mh will 
be increased by about 160 GeV because of the inclusion of the NLO corrections. 

• As illustrated in Fig lTTl the allowed region of Mh in model III will be shifted toward 
heavy mass end for a non-zero relative phase 9. 
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APPENDIX A: INPUT PARAMETERS AND RUNNING COUPLINGS 

In this appendix we list the input parameters used in our calculations. The values of 
most parameters are quoted from Ref. 0,01 directly. The masses m t , m b , m c are understood 
to be the pole mass, while rn q (fi) is the running q-quark mass in the modified minimal 
subtraction (MS) scheme at the renormalization scale /i. To first order in as, the running 
mass m q (fi) and the pole mass m q are related through 




(Al) 
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The masses, coupling constants and other input parameters are (H, 0] 



M Bd = 5.279GW, m b = 4.8 ± 0.2GeV, m c /m b = 0.29 ± 0.02; 

M w = 80A2GeV, m t = 174.3 ± B.lGeV, 

a s (M z ) = 0.119 ±0.004, a em = 1/137.056, 

Aj = -0.50GeV 2 , A 2 = 0.12GeV 2 

The ratio of CKM elements | V t *V t b/V c b\ 2 appeared in the decay rate B{B 
be expressed in terms of the Wolfenstein parameters as follows 



(A2) 
Xs^) can 



v; s v tb 



v, 



cb 



1 - A 2 (l - 2p) + A 4 (p 2 + f - A 2 ) + 0(X 6 ) 
0.971 ±0.010, 



where we have used A = 0.2196, A = 0.854, p = 0.22 ± 0.10 and fj = 0.35 ± 0.05 
error induced by 5f] = 0.05 is only 0.0001 for the ratio. 

For the semileptonic branching ratio Bsl = B(B — > X c eV e ) we use 

B SL = (10.64 ±0.23)% 



(A3) 
The 

(A4) 



as given in Ref. 0. 

Finally, for the running of a s (p), we use the two-loop formulae as given in Ref . jlol] : 



with 



a s (M z 
v(p) 



v(p) 



1 - 



A a s (M z ) ln(v(p)) 

P 47T v(p) 



a s {M z ) (M z 
Po — ~ m 



2tt \ p 

where /? — 23/3 and (3i = 116/3 for b quark decays. 



(A5) 



(A6) 



APPENDIX B: NLO WILSON COEFFICIENTS AT p = M w 

For the completeness, we list here the NLO functions at the matching scale pw ap- 
peared in Eqs. (j54l56l57J) . For more details see Ref.j3|. 
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In the SM, we have 
E (x) 



W 8 ,sm(x) 



T 7 ,sm(x) 
T 8 ,sm{x) 



x (x 2 + llx - 18) x 2 Ux 2 - 16a; + 15) , 2, 2 
12(x-l)3 + 6(^1? In^-gln^--, 

-16a; 4 - 122a; 3 + 80a; 2 - 8a; T . / 1 

Li 2 1 1 - 



X 



+ 



9(x- l) 4 
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3(x- l) 5 
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81(x- l) 5 
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+ 
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3 
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(Bl) 



(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 
(B9) 



where x = mf/M^. 

In the 2HDMs, the NLO functions describing the charged Higgs contributions are 
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M 7> yy(v) 
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FIG. 1: The box diagrams for the B° d -B° d mixing in the framework of 2HDM at the lowest order. 
More crossed diagrams are not shown. For those Feynman diagrams of NLO QCD corrections 
obtained by connecting any two quark lines with a gluon line, see the figures in Ref.jfJ]. 
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FIG. 2: Contour plot of \V t d\ vs the non-perturbative parameter fs d \J ' Bs d by using the well 
measured mass splitting AMs d = 0.502 ± O.OOTps" 1 . The shaded band corresponds to the limit 
Wtd\ = (7.9 ± 1.7) x 10~ 3 obtained from the measured AAf^ d . The width of the curve shows the 
effect of uncertainties of other quantities in Eq. ()ll|) . 
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FIG. 4: Plots of the mass splitting AMg d vs \Xu\ in model III for fB d \ BB d = 0.19 (the lower 



three curves) and fs d y Bb a = 0.23 (the upper three curves) and for Mh = 200 (solid curves) 
250 (short-dashed curves) and 300 GeV (dashed curves), respectively. 




FIG. 5: Contour plot of the mass splitting AMs d in the Xu — Mjj plane. The measured 
AMB d = 0.502 ± 0.007 is used. The A area is allowed by considering the error of AMs d only, 

but both A and B area will be allowed if the uncertainty of is d y Bs d is also taken into account. 
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FIG. 6: Plots of the mass splitting AMs d vs tan/3 in model II for fB d \ Bs d = 0.19 (the lower 



three curves) and fs d y Bs d = 0.23 (the upper three curves) and for Mh = 200 (solid curves) 
250 (short-dashed curves) and 300 GeV (dashed curves), respectively. For details, see text. 
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FIG. 7: Plots of M H dependence of C^ SM (M W ) (solid line) and the four functions Cj' YY (M w ) 
(the dashed and dot-dashed curves ) and C®' XY (Mw) (the short-dashed and dots curves). The 
decoupling behaviour of new physics contributions can be seen clearly. 
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FIG. 8: Contour plot of the branching ratio B(B — > X s j) in the |A#| — |A&&| plane, assuming 
M H = 250 GeV. 




FIG. 9: The M# dependence of branching ratio B{B — > -^7) for the Case A: (A^, A^) = (0.5, 1) 
and at NLO level. The band between two horizontal dot-dashed lines shows the data at la level. 
The dashed line and solid curve refers to the SM and Model III prediction, respectively. 



36 



6 




200 400 600 800 1000 

M H [ GeV] 

FIG. 10: The Mh dependence of branching ratio B(B — ► X s ^y) for the Case B: (Xtt,\b) = 
(0.5, 12) and at LO level. The band between two horizontal dot-dashed lines shows the data at 
2a level. The dashed line and solid curve refers to the SM and Model III prediction, respectively. 




FIG. 11: The Mh dependence of Wilson coefficient Cj (mi,) and its SM and model III parts 
for (Xu,Xbb) = (0.5,22) (i.e. the Case C).For details, see the text. 
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FIG. 12: The Mh dependence of the terms appeared in the curly brackets of Ea. (|69|) and the 
summation R in the SM and model III of Case C. For details, see the text. 
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FIG. 13: Plot of the branching ratio versus Mjj in model III of Case C. The band between 
two horizontal lines shows the data with 2a errors, the shaded band shows the SM prediction: 
B(B X s ~i) = (3.52 ±0.28) x 10~ 4 . The dashed and solid curve shows the LO and NLO model 
III prediction, respectively, (b) is a magnification of the light Higgs region of (a). 
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FIG. 14: Contour plot in |A&&| — Mh plane obtained by using the LO and NLO model III 
predictions and the measured decay rate at 2a level. The regions between two dashed curves 
and two solid curves are allowed when the LO and NLO theoretical predictions are employed. 
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FIG. 15: Plots of the branching ratio B{B — » X s j) vs the mass Mjj in model II for tan /3 = 4. 
For details see text. 
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FIG. 16: Contour plot in tan/3 — Mh plane obtained by using the NLO model II predictions 
and the measured decay rate at la (dashed curve) and 2a level (solid curve). The excluded 
region is below the corresponding curves. 
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FIG. 17: Plots of the branching ratio B(B — ► X s j) vs the mass Mu in model III for case C 
and for 9 = 0° (solid curve), 30° (dots curve), 60° (short-dashed curve) and 90° (dashed curve), 
respectively. The band between two horizontal dot-dashed lines shows the measured branching 
ratio within 2a errors. 
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FIG. 18: Plots of the branching ratio B(B — > X s ^) vs the phase 6 in model III for the case 
A (short-dashed curve ) and case C (solid curve), and assuming Mh = 250 GeV. The band 
between two horizontal dot-dashed lines shows the measured branching ratio within 2a errors. 
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